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K^ ■ Abstract. 

•y^ ' For a dimer with a non-degenerate orbital built from atomic wave functions 

of Gaussian shape we evaluate all the electron-phonon couplings derived from 
04 ■ the one-body and two-body electronic interactions, considering both the adi- 

0^ ' abatic and extreme non-adiabatic limit. Not only the values of the coupling 

^D ' parameters in the two limits, but also the expressions of the corresponding 

^D , terms in the Hamiltonian differ. Depending on the distance between the dimer 

^D . ions, some of the two-body couplings are comparable, or even larger than the 

one-body ones 

a 

^ : 1 The Model. 

n , 

>^ ' In a general two-site electron-phonon Hamiltonian H — H^i + Hph -\- Hei-ph the 

^ ' interacting term Hd-ph originates from developing H^i to first order in the ion 

displacements Ui (i = 1,2), where, in standard notation for a non-degenerate 
ij orbital: 

H \ Hei ^ e'^{nicr+n2a)+'^[t+X{ni^„+n2-cr)]{c\„C2a+II.c.)-{-U{ni^nii+n2'\n2i) 



+ {V-J/2)nin2~2J 



S^iS^ + ^{S+S2+H.c.) 



P{cl^cl^c^^C2^+H.c.). (1) 



We shall develop both the one-(e,t) and two-body ( A, t/, V, Jz = Jxy ~ P) 
electron interaction parameters, evaluated as in Rcf.l, by assuming a non- 
degenerate orbital described by Wannier functions built from atomic orbitals of 
Gaussian shape. We associate to each site a Gaussian atomic- like orbital 4>i{r — 
Ri), with the ions centered at the positions Rj = (±a/2 -|- Ui, 0, 0) (i = 1,2). 

By defining A EE (2/7r)^/'*r3/^ they read: 0j(r) = Nexp j-P^ (x ± a/2 - u;)^ + y 

Then the Wannier functions ^i, ^2 can be written as: 

*i(r) = A(S)0i(r) + B(S)02(r) *2(r) = B(S)0i(r) + A(S)02(r) 



2 



A{S) = [(1 + S)-^/^ + (1 - S)-^/^]/2 B{S) = [(1 + S)-^/^ - (1 - S)-^/^]/2 

(2) 
with S{u) = ((/)i|02) = exp[— r^(a + u)^/2], and m = U2 — ui. To make clear 
our method of calculation, it is convenient to explicitate, as an example, the 
one-body local electronic energy: 



e«= / vI/*(r,Ri,R2) 



-— V^ + Vi(r - Ri) + V2(r - R2) 



*i(r,Ri,R2)d^r 



..(3) 
where the potentials originating from the ion cores at the displaced positions 

Ri and R2 are: 

K, = y(r-Ri) = -e2z[(x±a/2-Ui)Vy^+z2]"'^^ (i = l,2) (4) 

with — e the electron charge, and +Ze the charge of the ion core. The local 
energy e (actually site-independent) can be decomposed into three terms, corre- 
sponding to the contributions from the kinetic operator (ey ) and from each one 
of the ionic potentials (ey , ey , respectively) . A similar decomposition holds for 
the hopping amplitude t. 

We need to distinguish between the adiabatic and the non-adiabatic limit 
in evaluating the electron-phonon interactions, because the integrals have dif- 
ferent kernels in the two cases. Indeed, in the adiabatic limit the displacements 
affect both the potentials V{r — Ri) and the Wannier functions ^j(r, Ri,R2) 
, expressing the requisite that the electronic charge distribution adjusts itself 
instantaneously at the position of the ions. We shall schematize the opposite 
situation, where the electrons are slower than the ions, as realized by the elec- 
tronic charge distribution staying centred around the undisplaced ion position, 
while the potentials are centred on the displaced ions. We shall call this the 
extreme anti-adiabatic limit. In the literature [gl the two limits are also named 
from, respectively, Frolich and Bloch. 



2 Coupling terms in the anti-adiabatic limit. 

In the anti-adiabatic limit Cy does not change, therefore no electron-phonon 
coupling originates from it. The coupling terms derived from the two-body in- 
teractions are also identically vanishing in this limit, because they involve the 
Wannier functions and the inter-electronic Coulomb potential which are both 
insensitive to the displacements of the ions. The only non- vanishing electron- 
phonon non-adiabatic interactions arise from the variation of the potential con- 
tributions to e and t. Let us now succinctly describe their evaluation. Full 
details are given in Ref.3. 

(i) , , 

The perturbation of Cy originates a term in the Hamiltonian connecting the 

local charge with the local deformation: X]cr(5o "-lo-'^^i + 5o "-20-^2) . It has 
a formal similarity with the Holstein coupling term[H, but it has a physically 
different origin, as HolsteinB considered electrons moving along a chain of fixed 
spacing with vibrating diatomic molecules at its nodes. On site 1, its explicit 



evaluations yields: 



^(1) ^ 2T^ 1^2 [Fo{2^V^) - S^] + AABS 



ftl^l-5 



(5) 



where Fo(x) — x ^/^Erf{x^^^). As, under site permutation, a —^ —a, ^q = 

(2) 

— 5q , as might have been anticipated by considering that, for equal charges 
ni — n2 and displacement amplitudes, with ei2 = (R2 — R'l)/^-, the energies 

gQ riiUi • ei2 and ^q n2U2 • ei2 on both sites coincide. Now symmetry requires 
Ui = — U2 (a constraint which does not hold for the original Holstein model[H) 

from which ^q — —g^ follows. The contribution to the Hamiltonian then 
reads: g^^^ EaC'^i.rWi - '^2<t'"2) 

The "crystal-field" coupling term g),i expresses the change in the energy 

Cy . To establish the form of this term, let us consider site 1, with charge ni. 
Its energy, after a displacement U2 of the ion on site 2, changes by an amount 
-E'^^ — gL 'ni U2 • ei2. This can be considered as the quantity measured by an 
observer sitting on ion 1 and watching the ion 2 moved by U2 . The equivalent 
measurement done by an observer on ion 2 watching the ion 1 displaced by Ui , 

yields E^'^' ~ g^r n2Ui • ei2. Assuming ni — n2,andui — ~~U2 one must have 
i?(^^ = £''^', implying g\.f = —g\.f ■ Therefore for the dimer as a whole one 
writes this term as g^? X]<j("'io-W2 — n2aUi). The explicit evaluation[pl for site 1 
yields: 



(12) 



= -2A (-\ J- [v4Fo(2a2r2) + ABSF^ia^V^ /2) - 4BS^ - AS^] . (6) 

From EqJq and the change of sign of a under site permutation, g^^^-* — — g^^^-* 
follows. 

The Su-Schrieffer-Heeger[H (SSH) interaction 7(^'^\ characterizing the SSH 

Hamiltonian Hssh = 7*-^^^ S<t("^1o-'^2ct + C2crCicr)(u2 — ui), is due to the modu- 
lation of the hopping amplitude t. To preserve the invariance of Hssh under 
site permutation, 7^^^^ = —7(21) ^ggg g g ^gfg 5 a^fj jy Indeed, its exphcit 
expression 1^ is: 

+ 4^1 (^^) {{A' + B^)So [So - Fo{a'ry2)] } . (7) 

Under site permutation a -^ —a and A ^>- B so that 7'^^) = —7(21) as expected. 
In conclusion, in the non-adiabatic limit the complete electron-phonon Hamil- 
tonian is given by: 

a (7 

+ 7^''^ E(^l-^2. + cLci J(^2 - ui) (8) 



3 Coupling terms in the adiabatic limit. 

In the explicit expression of the different electronic interactions in the adiabati- 
cally displaced state, u invariably enters in the combination a + u. So, to obtain 
the corresponding couplings, one can simply take the derivative with respect to 
a of the parameters in Eq^yJ as evaluated in Ref. 1 . 

There is some confusion in the literature about the correct form of the 
electron-phonon Hamiltonian obtained, in the adiabatic limit, from the vari- 
ation of the local energy e, therefore we shall devote some space to clarify this 
point. In this limit, Ui ^ Q enters both the charge distributions and the po- 
tentials. The " crystal field" interaction couples the charge on site i to the 
relative position of site j through the modification of both the kinetic and 
the potential contributions to e^*^. Let's place the origin of the a:-coordinate 
onto one of the displaced ions, at Ri, say. One has then to take into ac- 
count the relative displacement of the ions. In the adiabatic limit, there- 
fore, the overall e-derived electron-phonon coupling term in the Hamiltonian 

(12) v^ / N , (21) v^ / \ K (12) (21) 

is: ^ gi 2^^ nio-(u2 - "i) + gl ^ l_,„ "-2o-(wi - U2). As gi ' = -g\ we can 
write the total adiabatic contribution from local energy terms to the electron- 

phonon Hamiltonian as ii%^ — g\ X]<t(^2o- + ^1(t)(m2 — ""i)- A Hamiltonian 
of this form was used in the papers of Ref. 8, while those of Ref. 9 proposed 
Hamiltonians incompatible with our results. After including the SSH term, 
the complete one-body electron-phonon Hamiltonian in the adiabatic limit has 
therefore the form|[10|: 



Kp - gi''^ E("l- + "2-)("2 - ^1) + 7^''^ E(^l'^C2a + cLci J(W2 " ^l). (9) 
a a 

Coming now to the explicit expressions of the one-body coupling parameters, 

1 i -i r ■ (12) (12) , (12) J (-12) (12) , (12) -jjr 

let us write tor convemency gl ' = g\j ' + gy and 7'-^'''' = 7y + 7y ■ We 
obtain: 



(12) 



(12) 



[2 (1 - a^V^ - S^)\ 





du a 

du a 



(10) 



7v 



(12) 



2m 



aST^ 



(1-52)2 



[2{l-S^)-a^r^{l + S^ 



(12) 

Tv 




dAB , ^^±^g2 ^ AB^^ 



S 



du a 

djA^+B"^) A^ + B^ 
du a 



(l + a^T^) 



dAB AB 



du a 



i^o(a'rV2) 



Fo{2a^T^) 
(11) 



Notice that, as the partial derivatives are hnear in a, they change sign under 
site permutation. 

Also the two-body electronic interactions U,V,J{~ P),X of Eq. (1) give 
rise to electron-phonon couplings, all of the form: 



dY 



da 



H?-{-r-] Ficl,c^,){u, - u,) (*,j = 1,2) 



(12) 



where Y — U,V,X,J and F{Ci^,Cja) is the function of Fermi operators rep- 
resenting the two-body interaction whose amplitude is Y. From the results of 
Ref.l, their evaluation is trivial. Notice that, as U{a) = J(a) + e^Y j ^pH^ then 
dU/da = dJ/da = dP/da. We list below their explicit expressions: 



da 



i-aTS) 



S/a 



(1 + 352) 



(1-^2)2 

o r 



(1-52)3 



4a2r252 



1 + 252 + Fq (a^r^) - 2(1 + S^)Fo ( 
S^-Fo{a^T^)\ 



a2r2 



F„,^,-l 



S/a 



(1-52)2 



2(1 + 5^) 



^0 



a2r2 



(13) 



dU 
da 



-4aT^S^ 



+e^ 



(1-52)3 
52/a 



2 - 52 + 25"* + 52^0 (a2r2) - 452i^o 



a2r2 



2a^T^ 



(1-52)2 
+ S^-Fo {a^T^) + 4 



1 - 452 - Fo (a^T^ 



4Fn 



a2r2 



^o(^]-V5 



(14) 
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r 4ar252 1 




" 








;-; 


V^ 


[ (1-52)3] 






+e'^ 


1/a 




0? 


[(1-5 


2 


'J 



diV-J,/2) 
da 



3 - 52 - 85^ - (7 - 552)fo(a'r2) - 4(1 - 352)fo 



2a^T^S^ 



-1 - AS^ + Fo (a'^r^) +4:Fo 



F„i°^]--/S 



55 


v.( 


a2r2\l 
4 jj 


/a2r2\] 




I 4 )\ 




5 


>, 


(15) 


J 


, 





+ 252 + ^54 - (2 - hA ^o(a'r2) + 252 

Fig. 1 and 2 present the values of the one-body coupling constant for, respec- 
tively, the non- adiabatic, and the adiabatic limit, evaluated by assuming for 
the shape-controlling parameter of the Wannier functions F the typical R] value 

l.A~^. The most unexpected result concerns g^j- While usually neglected in 
the literature M on metallic systems, this coupling has been recognized as rele- 
vant to polar materials ||ll[. We find indeed that, when F = l.oA"^, g^-f is larger 
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Figure 1: Non-adiabatic coupling constants g^ ,g]^f ,7''^^-' (in cVA ^) versus 
the dimer length (in A), evaluated assuming F = l.oA^^. 
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Figure 2: Adiabatic coupling constants g): ,7*-^^-' (in eVA ^) versus the dimer 
length (in A), for F = 1.0 A^i. 
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Figure 3: Adiabatic coupling constants from two-body interactions dX/da, 
dU/da, d{V - J/2)/da (in cVA^i) versus the dimer length (in A), for T = 1.0 



than go for any a, and it becomes the largest parameter for a > 2.2A. For small 
a, the SSH coupling is the largest. 

(12) 

Fig. 2 for the adiabatic case shows that g^ is always larger than the SSH 
interaction 7*^^^% and particularly for large a there is an order of magnitude 
difference between them. The case of F = 2.A~^ (discussed in ref.3) shows 
that, the more localized are the orbitals, the more relevant is the role of g^/ in 
relation to the other admissible couplings. Fig. 3 shows the couplings derived 
from the two-body electronic interactions for the same parameters as Fig. 2. 
In general, their values are much smaller than those of (/e and 7'-^^', with the 
possible exception of d{V—J/2)/da. Indeed, that couphng arises from a physical 

mechanism not very different from the one originating gy, , i.e. the vibration 

of the charge on site j as felt by site i. Similarly to g^ also d{V — Jz/'i^/da 
decreases slowly with a, so that for large a those two are the only relevant 
couplings. Such interactions in the lattice have been recently discussed in ref.l2, 
while their effects in the optical spectra have been treated in ref.l3. 



4 Conclusions 

We have presented the analytical evaluation of the clcctron-phonon coupling 
parameters derived from both one- and two-body electronic interactions in a 
model of a dimer with a non-degenerate orbital built from atomic orbital of 
Gaussian shape. We have shown that the coupling terms in the adiabatic and 
the anti-adiabatic limits differ qualitatively. 

The evaluation of the coupling terms originating from the two-body elec- 
tronic interactions shows that at least the one generated by the Coulomb repul- 
sion between the charges on different sites, is comparable to, or even larger than, 
the couplings derived from the one-body interactions. The quantitative results 
for the coupling parameters, even if agreeing in order of magnitude with some 
estimates from experimental data|14 are obviously model-depending. However, 



their ratios should be more close to the reality. In particular, the obtained values 
of the coupling terms, when compared to the values of the electronic interactions 
resulting from the same Wannier functions [0 suggests that, for dimer lengths 
comparable to the lattice parameters in high temperature superconductors and 
colossal magnetoresistance materials, at most dU /da and dX/da can be safely 
dropped, while neglecting any of the other electron-phonon interactions is a 
questionable approximation. 
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